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Abstract 

A set of lower bounds on the continuum percolation threshold rj c of overlapping convex hy- 
perparticles of general nonspherical (anisotropic) shape with a specified orientational probability 
distribution in d-dimensional Euclidean space M. d have been derived [S. Torquato, J. Chem. Phys. 
136, 054106 (2012)]. The simplest of these lower bounds is given by rj c > v/v ex , where v ex is 
the d-dimensional exclusion volume associated with a hyperparticle and v is its d-dimensional vol- 
ume. In order to study the effect of dimensionality on the percolation threshold rj c of overlapping 
nonspherical convex hyperparticles with random orientations here, we obtain a scaling relation 
for 7/ c that is based on this lower bound and a conjecture that hyperspheres provide the highest 
threshold among all convex hyperparticle shapes for any d. This scaling relation also exploits 
the recently discovered principle that low-dimensional continuum percolation behavior encodes 
high-dimensional information. We derive an explicit formula for the exclusion volume v cx of a 
hyperparticle of arbitrary shape in terms of its d-dimensional volume v, surface area s and radius 
of mean curvature R (or, equivalently, mean width). These basic geometrical properties are com- 
puted for a wide variety of nonspherical hyperparticle shapes with random orientations across all 
dimensions, including, among other shapes, various polygons for d = 2, Platonic solids, sphero- 
cylinders, parallepipeds and zero-volume plates for d = 3 and their appropriate generalizations for 
d > 4. Using this information, we compute the lower bound and scaling relation for rj c for this 
comprehensive set of continuum percolation models across dimensions. We demonstrate that the 
scaling relation provides accurate upper-bound estimates of the threshold rj c across dimensions and 
becomes increasingly accurate as the space dimension increases. 

PACS numbers: 64.60.ah, 05.70.Jk, 71.30.+h, 82.70.Dd 
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I. INTRODUCTION 



Clustering and percolation behavior of many-particle systems is of relevance in a host of 
physical and chemical phenomena, including gelation and polymerization, structure of liq- 
uids and glasses, hopping in semiconductors, metal-insulator transition in condensed matter 
systems, nucleation, condensation of gases, chemical association, conduction in dispersions, 
aggregation of colloids, and flow in porous media Qfl. In a previous investigation fl, 
analytical estimates of the percolation threshold n c for continuum (off-lattice) percolation 
models of overlapping hyperparticles in ci-dimensional Euclidean space lR d were obtained for 
any d. An important consequence of this analysis, which we exploit in the present work, is 
that the large-ci percolation value for any hyperparticle shape is an important contribution 
to the corresponding low- dimensional percolation value. In other words, low-dimensional 
percolation properties encode high-dimensional information. The analysis was aided by a 
remarkable duality between the equilibrium hard-hypersphere (hypercube) fluid system and 
the continuum percolation model of overlapping hyperspheres (hypercubes), namely, 

P(r;r}) = -h(r;-rj), (1) 

where P(r; rf) is the pair connectedness function at some radial distance r, the quantity 

77 = pv (2) 

is a dimensionless or reduced density, p is the number density and v is the volume of a 
hyperparticle js]. Additionally, h(r; rf) is the total correlation function for the corresponding 
equilibrium hard-particle models at packing fraction r), fraction of space covered by the hard 



particles. The duality relation (P) relates a statistical property of a topological problem on 
the one hand to that of a geometrical problem on the other It is exact in the large-rf 



limit and is a relatively accurate approximation in low dimensions 



It was also shown in Ref. 
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a. 



that lower-order Pade approximants of the mean cluster 
number S of the form [n, 1] (where n = 0, 1 and 2) provide upper bounds on S and hence yield 
corresponding lower bounds on the threshold n c for ci-dimensional overlapping hyperspheres 



and overlapping oriented convex hyperparticles with central symmetry M (e.g., spheres, 
cubes, regular octahedra, and regular icosahedra for d = 3) [7| . The simplest of these lower 
bounds on rj c ([0, 1] Pade approximant) is given by 

1° * & ( 3 ) 



where it is to be noted that 2 d is the ratio of the exclusion volume v ex associated oriented 
convex hyperparticle (which includes the hypersphere) to the volume v of the hyperparticle. 
This trivially translates into a lower bound on the mean number of overlaps per sphere at 
the threshold M c : 

K EE 2% > 1. (4) 

The lower bound estimates as well as Percus-Yevick-like approximations were demonstrated 
[?| to become accurate in relatively low dimensions, improve in accuracy as d increases, and 
become asymptotically exact in the large-d limit, i.e., 

Vc d-^oo (5) 

and 

M c ~l, d^oo. (6) 

The aforementioned trends and asymptotic results were shown to hold for overlapping 
convex d- dimensional particles (or hyperparticles) of arbitrary shape and orientational dis- 
tribution when appropriately generalized 7|. For example, for identical overlapping hyper- 
particles of general anisotropic shape of volume v with specified probability distribution of 
orientations in d dimensions, the simplest lower bound on r\ c and on M c generalize as follows: 

'-'ex 

Mc EE Vc — > 1, (8) 

where 

w ex = / f(T,w)p(u)drdu), (9) 

f(r,u)) is still the exclusion-region indicator function, r and u) is the centroid position and 
orientation of one particle, respectively, with respect to a coordinate system at the centroid 
of the other particle with some fixed orientation, and p{u)) is the orientational probability 
density function. Moreover, in the high-dimensional limit for any convex hyperparticle 
shape, the following exact asymptotic results |13| immediately follow: 

v 

rj c ~ — , d — > oo (10) 

"^ex 

and 

JV C ~1, d^oo. (11) 



It was demonstrated in Ref. Q that the lower bound (J7J) on r] c improves in accuracy in 
any fixed dimension as the particle shape becomes more anisotropic. This property will be 
exploited later in Sec. ED to obtain an upper bound on rj c . 

We note that the bounds ([7]) and ([S]), and the associated asymptotic limits (fTUj) and 
( II ip also apply to overlapping zero-volume (d — l)-dimensional hyperplates with specified 
orientations in R d when the relevant parameters are appropriately generalized. While these 
"infinitesimally thin" hyperplates in M. d have vanishing volumes, their exclusion volumes are 
nonzero at number density p. However, an appropriately defined "effective volume" v e g must 
be chosen in order to define a reduced density rj of the form (J2J) for hyperplates. The choice 
for this effective hyperplate volume that we use henceforth is the volume of a (^-dimensional 
hypersphere of characteristic radius r, i.e., 

n d/2 

"* = F(57^T) r ' (12) 

where r is related to some characteristic length scale of the hyperplate and T(x) is the Euler- 
Gamma function. Thus, bounds (J7|) and dHJ still apply for hyperplates when v in expression 
for the reduced density rj is replaced by the effective volume (|T2|) . i.e., 

Vc > ^- (13) 

V e 

Application of the aforementioned continuum percolation results for nonspherical hyper- 
particles was only briefly explored in Ref. Q. The main objective of the present paper is to 
carry out such an investigation for a wide variety of nonspherical convex hyperparticles with 
random orientations across dimensions in order to ascertain the effect of dimensionality on 
the corresponding thresholds. This will be done by first ascertaining the exclusion volume 
of a wide variety of nonspherical hyperparticle shapes, and then applying the lower bound 
(J7|) on rj c and a scaling relation for r] c , the latter of which we obtain in the present paper. 

All of the results described above reveal that the topological problem of percolation in 
low dimensions becomes a purely geometrical problem in high dimensions. This behavior is 
explicitly manifested by the duality relation ([T]) and the fact the threshold r] c is determined 
solely by the exclusion volume associated with a hyperparticle in the high-c? limit [cf. (llOp ]. 
a purely geometrical characteristic. Exploiting the latter result, the principle that low- 
dimensional results encode high- dimensional information [/J, and a conjecture (postulated 
here) that hyperspheres provide the highest threshold among all convex hyperparticle shapes 



for any d, enables us to devise the aforementioned scaling relation for r\ c for randomly oriented 
hyperparticles of general nonspherical shape, which is also shown to bound rj c from above. 

The rest of the paper is organized as follows. In Sec. [TT1 we derive an explicit formula for 
the exclusion volume v ex associated with a nonspherical convex hyperparticle of arbitrary 
shape in terms of its basic geometrical properties. In Sees. II I II and IIV} we compute these 
geometrical properties for a wide variety of nonspherical hyperparticle shapes with random 
orientations in two and three dimensions, and for d > 4. A conjecture that systems of over- 
lapping hyperspheres provide the highest threshold among all convex hyperparticle shapes 
for any d is postulated in Sec. |VJ In Sec. IVI[ we derive a scaling relation for i] c and show 
that it is also an upper bound on r] c . In Sec. I VII} we compute the lower bound ((?]) and 
scaling relation for r] c for the aforementioned comprehensive set of systems of overlapping 
hyperparticles and compare the results to numerical data for rj c , including those that we 
obtain here for the Platonic solids. Finally, in Sec. IVIIIt we close with concluding remarks. 

II. EXCLUSION VOLUME FOR RANDOMLY ORIENTED OVERLAPPING 
NONSPHERICAL HYPERPARTICLES 

The exclusion volume t> ex associated with randomly oriented overlapping convex non- 
spherical hyperparticles plays a central role in this paper. In what follows, we derive an 
explicit formula for v ex in terms of the (^-dimensional volume v , surface area s and mean 
radius of mean curvature R (or mean width) of a hyperparticle of arbitrary shape. 

A. Radius of Mean Curvature (Mean Width) 

Consider any convex body K in d- dimensional Euclidean space M d . Let v and s denote 
the (i-dimensional volume and surface area of K, respectively. For example, in R 2 , ti and 
s are the area and perimeter, respectively, of the two-dimensional convex body or disk. In 
Euclidean space M. d , the parallel body associated with the convex body K at distance e 
is equal to the sum of K and a Euclidean ball of radius e. This operation preserves the 
convexity and compactness properties associated with K. The notion of a parallel body or 
set is well known in convex geometry and has been applied in the physical sciences, even 
for non-convex geometries; see, for example, the so-called dilation processes that has been 
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FIG. 1: (Color online) Illustration of the width w of a convex body in two dimensions. In M 2 , w(9) 
is function of the angle 6 between the normal n of the parallel planes (shown as red arrows) and 
the convex body. The mean with w is obtained by averaging w(9) over the angle 9, i.e., integrating 
w(6) over the angle 6 and dividing by 2tt. 

applied in the study of heterogeneous materials 

Il2|. 

The famous Steiner formula expresses the volume of the parallel body v e at distance e as 
a polynomial in e and in terms of geometrical characteristics of the convex body K, such as v 



and s. For example, for the cases d = 2 and d = 3, Steiner's formula 15j yields respectively 



v + se + Tie 



(14) 



and 



where 



v + se + 4irRe 2 + 



47T6 3 



(15) 



(16) 



is its radius of mean curvature and w its mean width, the latter of which is a lineal (one- 
dimensional) measure of its "size." Consider the convex body to be trapped entirely between 
two impenetrable parallel (d — l)-dimensional hyperplanes that are orthogonal to a unit 
vector n in W 1 . The "width" of a body w(n) in the direction n is the distance between 
the closest pair of such parallel hyperplanes, implying that the hyperplanes contact the 
boundaries of the body. The mean width w is the average of the width w(n) such that n 



is uniformly distributed over the unit sphere S 



d-l 



Figure [T] illustrates the concept of 



the width of a convex body in 
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The Steiner formulas (I14D and (ITS]) can be generalized for a convex body in any 
dimension d: 



v e = 

k=0 



where Wk are trivially related to the quermassintegrals or Minkowski Junctionals 15(. Note 
that Wo = Wi = s and = v s (d] e), where 

7T d/ V 

^ (d;a) = r(i + rf/2) (18) 

is the volume of a d- dimensional sphere (hypersphere) of radius a. Of particular, interest in 
this paper is the lineal characteristic in (ITTj) . i.e., the (d — l)th coefficient: 

W d - X = n(d)R } (19) 

where the radius mean of curvature R in M. d is related to the mean width w by f fT6|) . and 

n(i > = roT^y < 20 ' 

is the total solid angle contained in <i-dimensional sphere. 



B. Explicit Formula for the (/-Dimensional Exclusion Volume 

General formulas for the exclusion volume associated with identical randomly oriented 
nonspherical convex particles in two and three dimensions have long been known 
For example, for d — 2, it is known [l^ that 

s 2 

v cx = 2v + — (21) 

Z7T 

where s is the two-dimensional surface area or perimeter of the convex body. Note that the 
first term of 2v is trivial, accounting for the fact that two bodies are excluding one another 
and hence the exclusion volume must be at least 2v. At first glance, one might surmise that 
the appropriate three-dimensional expression for the exclusion volume involves a second 
term that is proportional to s 3 / 2 /(47r), where s is the surface area of the three-dimensional 
convex body. However, instead the correct three-dimensional formula 17| is 



S 



M — I I ( i- + i- )dS (23) 



where 

M= I , 

y 2 v#i i?a, 

is the integral of mean curvature (integrated over the entire surface of the convex body), 
and R± and R2 the principle radii of curvature. The latter quantity is related to the radius 
of mean curvature via the relation 

The second terms in the two- and three-dimensional formulas (j2ip and f )22|) . respectively, 
appear to be functionally distinct from one another. However, we will see shortly that by 
expressing both formulas in terms of R a unifying formula emerges that not only is valid for 
d = 2 and d — 3, but for arbitrary d. 

To our knowledge, a formula for the exclusion volume associated with identical randomly 
oriented d- dimensional convex nonspherical bodies K for any d has heretofore not been given 
explicitly. We find that this exclusion volume is given by 

2(2 d ~ 1 - 1) - 

v ex = 2v + -i — ! -sR. (25) 

We arrive at this general formula for any d by first recognizing that for the special case 
d = 2, the radius of mean curvature is trivially related to the two-dimensional surface area 
or perimeter, namely, s = 2itR, which enables us to rewrite relation ( 12TT) as follows: 



v ex = 2v + sR. (26) 

Similarly, using Eq. (|24|) . we can re-express relation (122"]) for d = 3 in terms of R, namely, 

v cx = 2v + 2sR. (27) 

We see that these two- and three-dimensional formulas differ by a coefficient multiplying the 
product sR, independent of the shape of the convex body. Recognizing that this coefficient 
depends only on dimension leads to the general d- dimensional formula for v ex associated 
with any convex body K that must have the form 

v ex = 2v + C(d)sR, (28) 

where C(d) is a d- dimensional coefficient that still must be determined. To determine C(d) 
as a function of d, we exploit the fact that formula ff28|) must apply to a d- dimensional sphere 



(hypersphere) for which v^/v, among all convex hyperparticles of nonzero volume, takes on 
its mnnmum value of 2" Q. For a hypersphere of radius a, the volume , is g,ve„ by 
R = a, and the surface area is given by 

s Jd; a) = dn ' .V " 1 . (29) 

1 ' ; r(i + |) 1 ; 

Using these last four geometrical-property relations in conjunction with formula (j28p yields 

2 d — 2 

C(d) = (30) 

Substitution of this expression into relation (128 j) gives Eq. (1251) for the exclusion volume. 

Using the general formula (1281) and the fact that v ex /v is minimized for hyperspheres 
(with v cx /v = 2 d ) among all convex hyperparticles of nonzero volume, enables us to obtain 
the following inequality involving s, R and v\ 

sR > dv, (31) 

where the equality holds for hyperspheres only. This bound is a type of isoperimetric in- 
equality and has been proved in the study of convex geometry using completely different 



techniques [18J. Note that for d = 2, inequality ( 13T|) yields the standard isoperimetric in- 
equality s > Attv, where we have used the fact that the surface or perimeter s = 2ttR. The 
c?-dimensional generalization of this standard isoperimetric inequality is 

s Q(d) d > dv^ , (32) 



where Q(d) is given by ( 120]) 18[. It is clear that this inequality is not the same as ( 13T|) . 

We see that in the case that a (d— l)-dimensional convex hyperplate in R d that possesses 
zero volume, i.e., v — 0, the exclusion- volume formula (12 5 j) yields 

^ex = 2 -^— r ^sR. (33) 
We shall consider examples of such convex bodies for d > 3 in the subsequent section. 



III. RADII OF MEAN CURVATURE AND EXCLUSION VOLUMES FOR VARI- 
OUS NONSPHERICAL PARTICLES IN TWO AND THREE SPACE DIMENSIONS 

Here we present results for the aforementioned geometrical characteristics of a variety 
of convex bodies that are randomly oriented in two and three dimensions possessing both 

10 



nonzero and zero volumes. In the former instance, it is convenient to consider the dimen- 
sionless exclusion volume v ex /v. Two- and three-dimensional exclusion volumes are known 



for a variety of convex nonspherical shapes that are randomly oriented lg,|l7j|. It necessarily 



follows that such exclusion volumes are larger than the corresponding values when all of the 
particles are aligns 
particle shape [7]. 



partic.es are ahgrred and ce„t ra .. y symmetric, in which case ,,>r = independent of the 



A. Two Dimensions 

Table H] provides volumes, surface areas, radii of mean curvature, and exclusion volumes 
for some two-dimensional convex particles. This includes the circle, certain regular polygons, 
rectangle, sphero cylinder and ellipse. The exclusion volumes for these particles are easily be 
obtained from relation (12 ip . Alternatively, the radius of mean curvature R can be obtained 
from the volume of the parallel body associated with the convex body at distance e that is 
of order e and use of relation (1191) . The corresponding exclusion volume can then be found 
employing relation ( )25|) for d = 2. 

In the "needle-like" limit for an ellipse {b/a — > oo), the general quantities R and v cx /v 
given in Table [J become 

- 2b v ex 8 b 

R ~ — , — ~ — --, 6/a-^oo. (34) 

The radius of mean curvature for a randomly oriented, zero- volume lines of length a in M 2 
is easily obtained from the rectangular case in Table U with a\ = and <22 = a, yielding 

R = - (line). (35) 

7T 

The corresponding exclusion volume is extracted from Table [I] in the rectangular case by 
first multiplying both sides of the relation for v ex / v by the volume v = ab, and then setting 
a\ = and 02 = a, which gives 

v ex = — (line). (36) 

7T 

Not surprisingly, R for a line of length a is identical to that of a needle-like ellipse of major 
axis 2b = a [cf. ([341]. 
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TABLE I: Volumes, surface areas, radii of mean curvature, and exclusion volumes for some convex 
particles K in M 2 . The radius of mean curvature and exclusion volume correspond to randomly 
oriented convex particles having nonzero volumes. 



K 



R 



Vex/v 



Circular disk, radius a 
Octagon, side a 
Hexagon, side a 
Pentagon, side a 

Square, side a 
Equilateral triangle, side a 
Rectangle, sides ai, (J2 
Spherocylinder, radius a 
cylindrical height h 
Ellipse, semiaxis a 
semiaxis b > a 



2(1 + \f2)o? 



-y/25+10y / 5a 2 
4 

a 2 

V3a 2 
4 

aia2 
7ra 2 + 2ah 
nab 



2na 
8a 
6a 
3a 

4a 
3a 

2(ai + 02) 
2(7ra + h) 



la 



3a 



5a 
2a. 



3a 
2tt 

ai+Q2 



§ + a 



4bE 



2 + 
2 

2 + 



16 



tt(1+V2) 

7T 

50 

V / 2^+lWf 



2 + 



2 + 



6^ 



2 + 



2(ai+a2) 2 



2 + 



2(7ra+fe) 2 
"a 



2 + -4£-£ 2 (,/l-2i ) 



B. Three Dimensions 

Table [III provides volumes, surface areas, radii of mean curvature, and exclusion volumes 
for some three-dimensional convex particles. This includes the sphere, Platonic solids, rect- 
angular parallelepiped, spherocylinder and spheroid. The radius of mean curvature R can be 
obtained from the volume of the parallel body associated with the convex body at distance 
e that is of order e 2 and use of relation ( JT9l) . The corresponding exclusion volume can then 
be found employing relation ( I25p for d = 3. 

It is easy to extract R and v ex /v from the spheroid results in Table ILT1 the "prolate-needle" 
(b/a —> 00) and "oblate-needle" (b/a — > 0) limits, which are respectively given by 

= b u ex 3tt b 

R~~, --, b/a -)> 00, 37 

2 v 4 a 

and 

R~-a, — ~- 7 , b/a->0. (38) 
4^46 
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TABLE II: Volumes, surface areas, radii of mean curvature, and exclusion volumes for some convex 
particles K in M 3 . The radius of mean curvature and exclusion volume correspond to randomly 
oriented convex particles having nonzero volumes. 



K 


V 


S 






Sphere, radius a 


3 


4-rro 2 


a 


8 


Icosahedron, side a 


i(3 + V5)a 3 


5v / 3a 2 


30 _ 1 , v/5 \ 
— cos 1 (^5)a 


90^ cos -1 f 

V 3 / 
^ 77(3+^) 


Dodecahedron, side a 


15+7^ a 


3^25 + lOv^a 2 


30 _i, 1 * 


/ " ^ 

90V25+10v / 5cos- 1 (4=) 
2_i_ v 5 


4 a 


^(15+7^) 


Cube, side a 


a 3 


6a 2 


3a 
4 


11 


Octahedron, side a 




2v / 3a 2 


3 —l/i \ 
2¥ cos (3)° 




Tetrahedron, side a 


t!« 3 


v^a 2 


Jjrcos-H— i)a 


2+ f \/i cos " 1( "^ } 


Rectangular 










parallclpiped 










sides a\, 0,2, 0,3 


010203 


JZi<j 2a i a 3 


1 v^3 

4 Z-/j=l a » 


/\ — '3 \ /x — *3 \ 

„ (£i<j«W»i)(£i=ia*) 
aia 2 a 3 


Cylinder, radius a 










height h 


na 2 h 


2na(a + fe) 


7ra+/i 
4 


n | (a+/i)(7ra+/i) 
^"i" 2ah 


Spherocy Under, radius a 










cylindrical height h 


ira 2 (4a+3h) 
3 


2?ra(2a + ft) 


ft+4a 
4 


„ . 3(2a+h)(4a+h) 
a(4a+3h) 


Prolate spheroid 










semiaxes a = c, b > a 










e 2 = 1 - (a/6) 2 


4-rra 2 6 
3 


27ra 2 (l + -^-sin-ie) 


f (1 + Fe tanh_1 e ) 


2+ |(1 + -^-sin-ie) 

1 2 v 1 ae / 

2 

X (1 + p-tanh _1 e) 


Oblate spheroid 










semiaxes a = c, 6 < a 










e 2 = 1 - (b/a) 2 


47ra 2 b 
3 


27ra 2 (l + J^tanh^e) 


+ ^sin- 1 ^) 


2 + Kl + ^tanh^e) 

X (l + ^ sin -i e ) 



Table HLT1 provides surface areas, radii of mean curvature, and exclusion volumes for some 
three-dimensional convex plates. This includes the circular, square, equilateral triangular, 
rectangular and elliptical plates. The radius of mean curvature R can be easily obtained 
from the volume of the parallel body associated with the convex body at distance e that is 
of order e 2 and use of relation (fl"9]) . The corresponding exclusion volume can then be found 
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using relation ( 125]) for d = 3. Alternatively, in the cases of circular, square and rectangular 
plates, R and v ex can be obtained from Table [Til as special cases of the cylinder (when h — 0) 
and rectangular parallelpiped ( = for the square plate and 03 = for the 

rectangular plate). To get t> ex in these cases, one must first multiply v cx /v in Table UTI by 
the volume v and then take the aforementioned limits. 

TABLE III: Surface areas, radius of mean curvature, and exclusion volumes for some convex plates 
K in R 3 possessing a volume v = 0. The radius of mean curvature and exclusion volume correspond 
to randomly oriented convex plates. 



K 


s 


R 




Circular plate, radius a 


2vra 2 


!« 


vr 2 a 3 


Square plate, side a 


2a 2 


\a 


2a 3 


Equilateral triangular plate, side a 


fa 2 


fa 


8 " 


Rectangular plate, sides 01, «2 


2ai<22 


\{ai + a 2 ) 


a\CL2{a\ + a 2 ) 


Elliptical plate, semiaxis a 








semiaxis b > a 


2-irab 


\bE{y/l - a/b ) 


2vra6 2 J B( v / l - a/b ) 



IV. RADII OF MEAN CURVATURE AND EXCLUSION VOLUMES FOR VARI- 
OUS HYPERPARTICLES 

We present here explicit formulas for the radii of mean curvature and exclusion volumes 
for a variety of hyperparticles in M. d with random orientations possessing both nonzero and 
zero volumes (hyperplates). 

A. Hypercube and Hyperrectangular Parallelpiped 

1. Hypercube 

The determination of the radius of mean curvature R for a <i-cube (hypercube) of edge 
(or side) length a is easily obtained from the d- dimensional generalization of Steiner's for- 
mula. Associated with the e neighborhood of an edge is a portion of a {d — l)-dimensional 
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hyper cylinder of height a and radius e, the volume of which is given by 

v cyl = v 8 (d-l;e)a= — - ^ ae \ (39) 

where v s (d;a) is the volume of a d-dimensional hypersphere of radius a explicitly given by 
f fTSj) . Hence, the volume associated with a hypercube edge is v cy y/2 d ~ 1 . Since a hypercube 
has a total of d2 d ~ 1 edges, the total volume associated with all of the edges is dv cy \. Therefore, 
equating this volume with the term Wd-it d ~ l in Steiner's formula ( TT71) and using (TI9~]) yields 

VL(d)Re d - 1 = dv cyh (40) 
and hence the radius of mean curvature for a hypercube of side length a is given by 

R = A ] l a . (41) 
For large d, we obtain the asymptotic result 

R ~ ^/J, d -> oo. (42) 

Not surprisingly, for large d, the longest diagonal, which grows like y/d, is the dominant 
contribution to R. 

Using Eqs. (I25j) and ( 14 ip and the fact that the surface area of a hypercube of side length 
a are given by s = 2d a d_1 and v = a d , respectively, we obtain the following expression for 
the exclusion volume: 

v ex 4(2«-i - 1) r(i + j) 



2 / 

In the \arge-d limit, this expression and ( 1421) yield the asymptotic result 

^L^ x f^2 d , d-*oc. (44) 



2. Hyperrectangular Parallelpiped 

Consider a <i-rectangular parallelpiped (hyperrectangular parallelpiped) with edges of side 
lengths ai, d2, • • • , ad- Using a similar analysis as for the hypercube, we find the radius of 
mean curvature is given by 

R= ^-±J) TV (45) 
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Employing Eqs. ( 125]) and ( l4"5j) and the fact that the surface area of a hyperrectangular 
parallelpiped of side lengths a\, 02, . . . , a<j is given by s = J2i<j<k<i--- 2> a i a j a kO>i • • • (where 
aidja k a[ • • • is the product of d — 1 different side lengths) and v = Yli=i a «' we obtain the 
following expression for the exclusion volume: 



2(2 



d-l 



^ 2a i a j a k ai ■■ ■ J ( ^ a» 

\i<j<k<l- 



1=1 



(46) 



« ' rf2 v^ r(^) nti 

Consider the special case in which there are d — 1 edges with side length a and the 
remaining edge with side length b. Then, formulas (145 p and (j4"6"I) simplify as follows: 



i2 



r(i + 



l)a + 6] 



^ex _ 4(2 d - 1 - 1) r(i + 
v d 2 ^ r(4± 



(d 



a 



(d 



(47) 



(48) 



We shall subsequently make use of these relations in deriving certain hyperplate results. 



B. Hyperspherocylinder 

Consider a d-spherocylinder (hyperspherocylinder) of height h and radius a, which can 
be decomposed into hypercylinder of height h and two hemispherical caps of radius a (i.e., 
hypersphere of radius a). The volume of order t d ~ x associated with the e neighborhood of a 
hyperspherocylinder is 

n (d-l)/2 e d-l h n d/2 de d-l a 

f(Si) + r(i + d/2)' 

which when equated with the term Wd-i^ 1 in Steiner's formula ( 1T71) and using ( 1T91) gives 

n(d)Re d - 1 = ^— - + — — , 49 

v ; r(^±i) r(i + d/2)' v ; 

and hence the radius of mean curvature for a hyperspherocylinder is given by 

R = 1. -g^r /i + a. 50 

dy?r(^±i) 1 ; 

For large d, we find the asymptotic result 



R ^ \l 7T-, h + a, d-^oo. (51) 
w 27rd 
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We see that if h grows mores slowly with dimension than \fd, 



R ~ a, 



d — > oo. 



(52) 



The surface area and volume of a spherocylinder are given respectively by 



s s (d — 1; a)h + s s (d; a) 



(53) 



and 

v = v s (d — l;a)h + v s (d;a), (54) 

where s s (<i;a) and v s (d;a) are the surface area of a <i-dimensional hypersphere of radius a 
given by ff29l) and ffl8|) . respectively. These relations combined with Eqs. (125]) and fl50|) yield 
the exclusion volume to be given by 



ex 
1> 



2 + 



2(2 



d-l 



d 



s s (d — 1; a)h + s s (c?; a) 
v s (d — 1; a)h + f s (rf; a) 



r(i + f)/i 



+ a 



(55) 



In the large-<i limit, this relation yields the asymptotic result 



h 



1 



a V 2nd 



1 



(56) 



where we have used (151j) . We observe that in order for the contribution to the the scaled 
exclusion volume from the cylindrical portion, first term in (|56|) . not to vanish as d becomes 
large, the aspect ratio h/a must grow as fast as y/d or faster. 



C. Hyperoctahedron and Hypertetrahedron 

A regular d-crosspolytope or regular hyperoctahedron (also known as the orthoplex) is 
the d- dimensional generalization of the three-dimensional regular octahedron. For d — 1 
and d = 2, the regular hyperoctahedron is a line and square, respectively. A regular d- 
simplex or hypertetrahedron is the ci-dimensional generalization of the three-dimensional 
regular tetrahedron. For d = 1 and d = 2, the regular hypertetrahedron is a line and 
equilateral triangle, respectively. For d > 5, there are only three types of convex regular 



polytopes: the regular hypercube, hyperoctahedron and hypertetrahedron 19|. For a regular 



hyperoctahedron and regular hypertetrahedron, the radius of mean curvature is not as easy 



to derive from Steiner's formula, but integral formulas for it have been derived 



20 



-221. 
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1. Regular Hyperoctahedron (d-Crosspolytope) 



For a regular hyperoctahedron with unit side length, the radius of mean curvature is 



given by 



where 



21 



22 1 



- y/2d(d-l) r £) 

R = y - i-^l i (d), 57 



poo 

Io(d) = / exp(-2x 2 ) erf(x) d - 2 dx. (58) 
Jo 

Note that for d = 2 and d = 3, we respectively have 



Io(2) = ^ (59) 



' o(3) = vi t *' 1_ ' ( 7! ) ' (60) 



and 



Substitution of these formulas into fl57|) recovers the results for R listed in Tables fl] and |TT] 
for the square and octahedron, respectively. For large d, it is straightforward to show that 

t tj\ ^Vlnrf 

W) 2^~' ^ 6 ^ 

which implies 



K-J^f, d-+oo. (62) 



The surface area and volume of a regular hyperoctahedron are given respectively by 

S ~ (d-l)! 

and 



(63) 



" = ~5P < 64) 

These relations combined with Eqs. (1251) and (1571) yield the exclusion volume to be given by 

v cx _ Ad^jd-i) r(j)(2^-i) 

V -2 + Io(d). (65) 

In the large-<i limit, this relation and (|62|) give the asymptotic result 

^/21n(d) 2 d , ci^oo. (66) 



v, 
v 



ex 
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2. Regular Hypertetrahedron (d-Simplex) 



For a regular hypertetrahedron with unit side length, the radius of mean curvature is 



given by 



where 



20 



22] 



R 



d(d+i) r(f) 



hid) 



exp(— 2x 2 



1 + erf(x) 



d-l 



dx. 



Note that for d = 2 and d = 3, we respectively have 



and 



/t(3) 



7r(2) 



1 



7T 



2^ 



cos 



1 



(67) 



(68) 



(69) 



(70) 



2\/27r v 3 

Substitution of these formulas into fIBTj) recovers the results for R listed in Tables [T] and [XX] 
for the equilateral triangle and tetrahedron, respectively. For large d, it is simple to show 
that 

27ry / in~d 

ir(d) j 2 — , (71) 

which implies 

/in A 

(72) 



i2 



d 



d — >■ oo. 



We see that in the large-ci limit, the radii of mean curvature for a regular hyperoctahedron 
and regular hypertetrahedron become identical. 

The surface area and volume of a regular hypertetrahedron are given respectively by 



and 



2(d-i)/2(d_ l)! 



(73) 



(74) 



2 d / 2 rf! ' 

These relations combined with Eqs. ( |25l) and (167}) yield the exclusion volume to be given by 

d^ 2 {d + \f' 2 r(f)(2 rf - 1 -i) 



'-'ex 
V 



2 + 



Trr(^i) 



-/r(d). 



In the large-d limit, this relation and (172"]) yield the asymptotic result 

— ~ ^dlnfd) 2 d , d^oo. 

v 
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(75) 



(76) 




FIG. 2: (Color online) The dimensionless exclusion volume v ex /v versus dimension d for the three 
convex regular polytopes: hypercube, hyperoctahedron and hypertetrahedron. 

Among the three possible convex regular polytopes in high dimensions (hypercube, hy- 
peroctahedron and hypertetrahedron), the hypertetrahedron has the largest dimensionless 
exclusion volume v ex /v for any fixed dimension d > 3, while the hypercube has the smallest 
value. Figure [2] shows v cx /v versus d for 3 < d < 11 for the three convex regular polytopes. 



D. Hyperplates 

1. (d- I) -cube in R d 

Consider a hyperplate in rf-dimensional space that is a (d — l)-cube with side length 
a. The radius of mean curvature and exclusion volume of such a hyperplate are trivially 
obtained from special cases of relations (j47p and (jUJ). That is, setting b = in (|47p yields 

(d-i)r(i + j) 

The corresponding exclusion volume is extracted from (|48|) by first multiplying both sides 
of the equation by the volume v = a d ~ 1 b, and then setting 6 = 0, yielding 



4(^-i)(2 d - 1 -i)r(i + f) 
rfV^TW) 



(78) 
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which, not surprisingly, is smaller than the exclusion volume ( fl3j) for a d-cube in M. d at fixed 
d. Not surprisingly, in the high-ci limit, formula (J77J for the radius of mean curvature tends 
to the same asymptotic expression (I42p for a G?-cube in M. d . This is not the case for the 
high-G? limit of (1751) for the exclusion volume, which is given by 

v ex ~ \hr 2d ad i d^oo, (79) 

V GOT 

and should be compared to (jSJ), which grows faster as d increases. 

Note that in the case d = 3, the formulas (1771) and (1781) give the same results reported 
in Table [Til for the instance of a square plate in IR 3 . For d = 2, they also yield the relations 
(135j) and (1361) for the case of a randomly oriented line in IR 2 . 



(d - l)-sphere in R d 

Consider a hyperplate in d- dimensional space that is a (d — l)-sphere of radius a. For 
such a hyperplate, R can be obtained from the corresponding relation (1771) for a cubical 
hyperplate by multiplying the latter by the factor 

yjr£(j) 

which is the ratio of radius of mean curvature R = a of a (g? — l)-sphere in ]R d_1 to that 
of the (d — l)-cube of side length a in the latter of which is obtained by using ( )4T|) . 

Thus, the radius of mean of curvature for a (d — l)-sphere of radius in IR d is given by 



- (d-i)r(|) 2 . . 

R = i >\ 2 > a . 80 

To illustrate the fact that this simple mapping is equivalent to the explicit calculation of R 
for a (d — l)-sphere, we carry out such a computation in Appendix A for a 3-sphere in M 4 . 
Using the fact that the surface area of a (d — l)-sphere in M. d is s = 2v s (d — 1; a), where 
v s (d; a) is given by (}T8~|) . we obtain the d- dimensional exclusion volume to be 

_ 2(ci-l)(2^-l) 7 r(^)/ 2 r(f) 2 2^/ 2 

dr(^ r(i + d/2) ' 1 j 

where the upper bound is the exclusion volume of a rf-sphere in M. d Not surprisingly, in 
the high-rf limit, formula ( IHUl) for the radius of mean curvature tends to the radius a. In this 
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asymptotic limit, relation ( )8T|) becomes 



23/2 f^m d/2 2 , ai< i f^mf 2 2 „ oi (82) 



ti v/7r V d J " " y/dn V d 
Observe that in the case d = 3, the formulas flSUj) and (jSTj) give the same results reported 
in Table HT1 for the instance of a circular plate in K 3 . For d = 2, they also yield the relations 
( 135|) and ( )36|) for the case of a randomly oriented line in M. 2 . 

V. CONJECTURE FOR THE MAXIMUM-THRESHOLD HYPERPARTICLE 
SHAPE AMONG ALL CONVEX BODIES IN R d 

Recall that the dimensionless exclusion volume v ex /v, among all convex bodies in M. d with 
a nonzero d- dimensional volume, is minimized for hyperspheres [see Sec. Ill B] and its thresh- 
old rj c exactly tends to v/v ex = 2~ d in the high-dimensional limit [cf. (J5J)]. These properties 
together with the principle that \ow-d percolation properties encode high-rf information 
leads us to the following conjecture: 

Conjecture: The percolation threshold n c among all systems of overlapping randomly 
oriented convex hyperparticles in M. d having nonzero volume is maximized by that for hy- 
perspheres, i.e., 

(Vc)s > Vc (83) 

where (rj c )s JS the threshold of overlapping hyperspheres. 

We note that similar reasoning leading to the aforementioned conjecture f )83|) for convex 
hyperparticles of nonzero volume also suggests that the dimensionless exclusion volume 
Vex/veE associated with a convex (d — l)-dimensional hyperplate in M. d is minimized by the 
(d — l)-dimensional hypersphere, which consequently would have the highest percolation 
threshold among all convex hyperplates. Recall that v e g the effective <i-dimensional volume 
associated with the zero- volume hyperplate of interest, defined by (Tl2l . 



VI. ACCURATE SCALING RELATION FOR THE PERCOLATION THRESH- 
OLD OF OVERLAPPING CONVEX HYPERPARTICLES 

Recall that lower bound ([7]) on the threshold n c of overlapping convex hyperparticles of 
general shapes with nonzero volumes and a specified orientational distribution is determined 
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by the dimensionless exclusion volume v cx /v, which also provides the exact high-c? asymptotic 
limit [cf. f fTUj) ]. A special case of this lower bound is the inequality (13 which is valid for 
hyperspheres as well as aligned centrally symmetric particles for which v ex /v = 2 d . 

Guided by the aforementioned high- dimensional behavior of rj c , conjecture ( 153"j) for hyper- 
spheres, and the functional form of the lower bounds ([3]) and (J7J), we propose the following 
scaling law for the threshold r] c of overlapping nonspherical convex hyperparticles of arbitrary 
shape and orientational distribution having nonzero volumes for any dimension d: 



where {f] c )s an d {v ex /v)s = 2 d are the percolation threshold and dimensionless exclusion vol- 
ume, respectively, for hyperspheres in dimension d. Thus, given the percolation threshold 
for the reference system of hyperspheres {rj c )s and the dimensionless exclusion volume v ex /v 
for some general nonspherical convex hyperparticle shape with a specified orientational dis- 
tribution, one can estimate threshold of such a system of overlapping hyperparticles across 
dimensions via the scaling relation ( I84p . Note that this scaling relation becomes exact in the 
high-<i limit. In the next section, we will show that this scaling relation provides reasonably 
accurate estimates of 1] C in low dimensions and, hence, must become increasingly accurate 
as d becomes large. This is yet another manifestation of the principle that low- dimensional 
percolation properties encode high- dimensional information [7]. 

It is noteworthy that the scaling relation (184")) not only provides a good approximation for 
the threshold rj c of overlapping nonspherical hyperparticles, it yields an upper bound on r) c , 
and hence a relatively tight one. Explicitly, we have the following d- dimensional inequality: 



This bounding property is a consequence of an observation made in Ref. \jX namely, that, 
at fixed d, the lower bound ([7]) on (r] c ) $ for overlapping hyperspheres converges more slowly 
to its exact asymptotic value of 2~ d than does the lower bound ([7]) for nonspherical hyper- 
particles with a specified orientational distribution. It then immediately follows from the 
functional form of ( I84"|) that it bounds f] c for nonspherical hyperparticles from above and 
converges to the exact asymptotic value of v/v ex in the large-c? limit [cf. f llOp ], 




(84) 




(85) 
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For a zero-volume convex (d— l)-dimensional hyperplate in M. d , it is more appropriate to 
choose the reference system to be (d — l)-dimensional hyperspheres of characteristic radius 
r, yielding the scaling relation for the threshold for some other hyperplate shape to be 

^ c ~ ( ~ ) ( ~ ) (^shp 

\V cS J SHP \V CX J 

= 2 d ( Vc ) SHP . (86) 

Here (r] c )sHP is the percolation threshold for a (d — l)-dimensional hypersphere and v e s is 
its effective volume, as defined by relation ( fl2l) . 

We will test these scaling laws and bounds for a variety of hyperparticle shapes across 
various dimensions using the exclusion volumes presented in Sees. IHII and HVl and comparing 
these results to available numerically-determined threshold estimates, including those values 
that we have obtained below for the Platonic solids. 

VII. APPLICATION OF THE SCALING RELATION AND BOUNDS 

Here, we apply the lower bound (JTj) and scaling relation (1841) to estimate r\ c for a wide 
class of randomly oriented overlapping nonspherical hyperparticles in IR d , including, among 
other shapes, various polygons for d — 2, Platonic solids, spherocylinders, parallepipeds and 
zero-volume plates for d = 3 and their appropriate generalizations for d > 4, when possible. 
We test these estimates against available numerical results for rj c in M 2 and M 3 , including 



those obtained here for Platonic solids using the rescaled-particle method 231 ] . 



Two Dimensions 



Table IIVI gives the percolation threshold i] c estimated using the scaling relation (jSi 
for a variety of randomly oriented overlapping particles in M 2 . In the cases of squares 
and ellipses with different asp ect ratios, the associated threshold value rf c obtained from 



numerical simulations 



2J28I 



-28] are also given in the table. Clearly, for most of randomly 
oriented nonspherical particles, the scaling relation not only bounds r] c from above, consistent 
with inequality ( 1851) . the bounds are relatively tight. Since the scaling relation already 
provides a good estimate for d = 2, it will increasingly become accurate as d increases for 
reasons discussed in Sec. ED We will see below that this is indeed the case for d = 3. 
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TABLE IV: Percolation threshold r\ c of certain overlapping convex particles K with random orien 
tations in R 2 , estimated using Eq. (|84p and the associated threshold values rj* for squares 



and ellipses 



circles 



24 



25] 



26-28] obtained from previous numerical simulations. The threshold of overlapping 



.l Jis also given, which is employed to compute the estimates. 



K 


V*c 


Vc 


Circle 


1.1281 




Square 


0.9822 


0.9925 


Ellipse b = 2a 


0.76 


1.107 


Ellipse b = 5a 


0.607 


0.7174 


Ellipse b = 10a 


0.358 


0.4402 


Ellipse b = 100a 


0.0426 


0.0543 


Octagon 




1.0980 


Hexagon 




1.0730 


Pentagon 




1.0463 


Equilateral Triangle 




0.8501 


Rectangle a<i = 2a\ 




0.9275 


Spherocylinder h = 2a 




1.0357 



B. Three Dimensions 



1. Nonzero- Volume Convex Bodies 



The percolation thresholds for randomly oriented overlapping Platonic solids (i.e., tetra- 
hedron, octahedron, dodecahedron, icosahedron and cube) are estimated using a recently de- 
veloped rescaled-particle method. Implementation of this algorithm for the Platonic solids is 
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the rescaled-particle method was applied 



described in Appendix B of this paper. In Ref. 
to obtain accurate estimates of rj c for overlapping hyperspheres and oriented hypercubes. 

Table |V] lists the numerical estimates of t] c , the associated lower-bound values t]l given 
by (|7j), and the numerical values of the dimensionless exclusion volumes for these polyhedra 
v cx/v (whose analytical expressions are given in Table HT|) . Observe that the tetrahedron and 
icosahedron possess the smallest and largest percolation thresholds, respectively, among 
the Platonic solids, which is confirmed in Table |V] As discussed above, this implies that 
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TABLE V: Numerically estimated percolation threshold rj c of randomly oriented overlapping Pla- 
tonic solids and the associated lower bound value tjl given by ([7]). The numerical values of the 
dimensionless exclusion volumes of these polyhedra v ex /v, whose analytical expressions are given 
in Table HU are also included 



K 




VL 


Vc 


Tetrahedron 


15.40743 


0.06493 


0.1701±0.0007 


Cube 


11 


0.09090 


0.2443±0.0005 


Octahedron 


10.63797 


0.09398 


0.2514±0.0006 


Dodecahedron 


9.12101 


0.1096 


0.2949±0.0005 


Icosahedron 


8.91526 


0.1126 


0.3030±0.0005 



the the tetrahedron and icosahedron possess the largest and smallest dimensionless exclu- 
sion volumes, respectively, among the Platonic solids. We note that our estimate of r) c 
for randomly oriented overlapping cubes is consistent with that reported in Ref. [24J, i.e., 
r] c = 0.2444 ± 0.0003, which verifies the accuracy of our approach for estimating rj c . The 
number within the parentheses represents the error in the last digit. These results will be 
employed below to verify the accuracy of the general scaling relation ( 1841) for estimating t] c 
of nonspherical convex hyperparticles across dimensions, since the scaling relation becomes 
increasingly accurate as d increases, as discussed in Sec. I VI I 

Estimates of the percolation threshold rj c using the scaling relation (l84j) for a variety of 
randomly oriented overlapping particles in M 3 are presented in Table IVII These predictions 
are compared to numerically determined threshold values. We see that for most of the ran- 
domly oriented nonspherical particles, the scaling relation ( 1841) . which is also an upper bound 
[cf. (1851)]. only slightly overestimates t] c compared with the numerical results. For oblate 
spheroids and prolate spheroid of aspect ratio 2, the scaling relation slightly underestimates 
r] c compared to the numerical values, implying that the numerical results are overestimates 
of the actual threshold for these systems, i.e., the numerical estimates in these instances 
violate the upper bound ( 1851) . These three-dimensional results show that scaling relation 
(I84p for T] c improves in accuracy over the two-dimensional results reported in Table IIVI 
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TABLE VI: Percolation threshold r/ c of certain overlapping convex particles K with random orien- 
tations in M 3 , estimated using Eq. (|84p and the associated threshold values r\* c for regular polyhedra 
(obtained from our numerical simulations) and spheroids 



23, 



28y . The threshold of overlapping spheres 



311 ] is also given, which is employed to compute the estimates. 



K 


* 
1 L, 


Vc 


Sphere 




0.3418 




Tetrahedron 




0.1701 


0.1774 


Icosahedron 




0.3030 


0.3079 


Decahedron 




0.2949 


0.2998 


Octahedron 




0.2514 


0.2578 


Cube 




0.2443 


0.2485 


Oblate spheroid a = c - 


= 1006 


0.01255 


0.01154 


Oblate spheroid a = c 


= 106 


0.1118 


0.104 


Oblate spheroid a = c 


= 26 


0.3050 


0.3022 


Prolate spheroid a = c 


= 6/2 


0.3035 


0.3022 


Prolate spheroid a = c 


= 6/10 


0.09105 


0.104 


Prolate spheroid a = c = 


= 6/100 


0.006973 


0.01154 


Parallelpiped 02 = 03 - 


= 2ai 




0.2278 


Cylinder h = 2a 






0.4669 


Spherocylinder h = 


2a 




0.2972 



2. Zero- Volume Convex Plates 



Table IVHI lists the percolation threshold r] c obtained from the scaling relation ( 1841) and 
numerical simulations for a variety of randomly oriented overlapping zero-volume two- 
dimensional plates in R 3 . As per relation (IT2"|) . the effective volume of these plates is chosen 
to be the volume of a sphere whose radius r is a characteristic length scale of the associated 
plate. In particular, for circular disks, r is chosen to be radius of the disk. For noncircular 
plates, r is the radius of a disk whose area is the same as that of the noncircular plate of 
interest [32] . Importantly, the reference system for the scaling relation is randomly oriented 
overlapping circular disks, which is also given in the table. Again, the scaling relation ((Hi 
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only slightly overestimates rj c in the case of the triangular and elliptical plates. For square 
plates, the scaling relation slightly underestimates r) c , suggesting that the numerical results 
overestimate the actual percolation threshold for this system. 



TABLE VII: Percolation threshold rj c of certain overlapping convex plates K with random orienta- 
tions in M 3 , estimated using Eq. (|84p . The associated threshold values rj* obtained from previous 



numerical simulations 



321 ] are also given if available. The threshold of randomly oriented over- 



lapping circular disks is also given, which is employed to compute the estimates. The numerical 
values of rf c listed here are given up to the number of significant figures reported in Ref. [32. 



K 


Vc 


Vc 


Circular disk 


0.9614 




Square plate 


0.8647 


0.8520 


Triangular plate 


0.7295 


0.7475 


Elliptical plate b = 3a 


0.735 


0.7469 


Rectangular plate 02 = 2eq 




1.0987 



C. Dimensions Four Through Eleven 



1. Nonzero- Volume Convex Bodies 



In Table fVllTl we provide the percolation threshold r\ c of (^-dimensional overlapping hyper- 
spheres and oriented hypercubes for 4 < d < 11 obtained previously via the rescaled-particle 
method [23]. As expected, the threshold values for the two systems approach each other as 
d increases and become identical in the limit d — > 00, i.e., r] c — > l/2 d [7]. The results for 
hyperspheres will be used to estimate rj c for the hyperparticles discussed in Sec. IV. 

To our knowledge, no numerical estimates of r] c are available for hyperparticles for d > 4. 
Simulations are difficult to carry out due to the increasing complexity in accurately detecting 
particle overlap as d increases. Given the accuracy of the scaling relation (1841 . verified 
by numerical results in M 2 and R 3 , we now apply it to estimate r] c for the d- dimensional 
randomly oriented overlapping hyperparticles discussed in Sec. IV, where it should provide 
even better predictions. Table IIXI summarizes threshold values for various hyperparticles 
for 4 < d < 11 using fl84"|) ) and the thresholds for hyperspheres listed in Table IVIIII As 
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TABLE VIII: Numerical estimates of the percolation threshold r/ c of d-dimensional overlapping 
hyperspheres and oriented hypercubes for 4 < d < 11 reported in Ref. 123 . The threshold values 



for the two systems, given up to the number of significant figures reported in Ref. l23j, approach 
each other as d increases and become identical in the asymptotic limit d — > oo, i.e., r/ c — > l/2 d 



Dimension 


Hypersphere 


Hypercube 


d=4 


0.1304 


0.1201 


d=5 


0.05443 


0.05024 


d=6 


0.02339 


0.02104 


d=7 


0.01051 


0.01004 


d=8 


0.004904 


0.004498 


d=9 


0.002353 


0.002166 


d=10 


0.001138 


0.001058 


d=ll 


0.0005530 


0.0005160 



TABLE IX: Percolation threshold rj c of certain d-dimensional randomly overlapping hyperparticles 
estimated using Eq. (|84p for 4 < d < 11, including hypercubes (HC), hyperrectangular paral- 
lelpiped (HRP) of aspect ratio 2 (i.e., a\ = 2a and Gtj = a for i = 2, . . . , d), hyperspherocylinder 
(HSC) of aspect ratio 2 (i.e., h = 2a), hyperoctahedra (HO) and hypertetrahedra (HT). 



Dimension 


HC 


HRP 


HSC 


HO 


HT 


d 


= 4 


8.097 x 10" 


-2 


7.452 x 10" 


-2 


1.109 x 10" 


-l 


6.009 x 10" 


-2 


3.471 x 10" 


-2 


d 


= 5 


2.990 x 10~ 


-2 


2.775 x 10" 


-2 


4.599 x 10" 


-2 


1.724 x 10" 


-2 


8.808 x 10" 


-3 


d 


= 6 


1.167 x 10' 


-2 


1.092 x 10" 


-2 


1.975 x 10" 


-2 


5.560 x 10" 


-3 


2.580 x 10" 


-3 


d 


= 7 


4.846 x 10" 


-3 


4.568 x 10" 


-3 


8.899 x 10" 


-3 


1.986 x 10" 


-3 


8.518 x 10" 


-4 


d 


= 8 
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expected, 7] c for randomly oriented hypercubes is lower than that for oriented hypercubes in 
the corresponding dimensions. 

The threshold rj c versus d for the systems listed in Table HXl is plotted Fig. El As expected, 
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FIG. 3: (Color online) Percolation threshold r\ c versus dimension d as obtained from scaling relation 
(I84j) for randomly oriented overlapping nonspherical hyperparticles, including hypercubes (HC), 
hyperrectangular parallelpipeds (HRP) of aspect ratio 2 (i.e., a\ = 2a and di = a for i = 2, . . . , d), 
hyperspherocylinder (HSC) of aspect ratio 2 (i.e., h = 2a), hyperoctahedra (HO) and hypertetra- 
hedra (HT). Values of r\ c for hyperspheres (HS) 23J are also shown for purposes of comparison. 

threshold values of nonspherical hyperparticles are always below that of hyperspheres for 
any fixed d. For hyp ersphero cylinders of aspect ratio 2 (i.e., h = 2a), the associated r] c is 
only slightly below the hypersphere value in the corresponding dimensions. Similarly, for 
hyperrectangular parallelpipeds of aspect ratio 2 (i.e., ai = 2a and = a for i — 2, . . . , d), 
T] c is only slightly below the hypercube value in the corresponding dimensions. 

To obtain excellent estimates of r] c for nonspherical hyperparticles for d > 12, one can use 
the scaling relation (1&4"1) together with highly accurate analytical expressions for hyperspheres 



in high dimensions given in Ref. 3; see, e.g., Eqs. (86) and (119) in that paper. 



2. Zero-Volume Convex Hyperplates 

The left panel of Fig. H] shows the dimensionless exclusion volume v ex /v e s versus d for 
spherical and cubical hyperplates. Recall that the effective volume v e s of a hyperplate [i.e., 
a (d — l)-dimensional object in M. d ] is chosen to be the volume of a d-sphere with radius 
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FIG. 4: (Color online) Left panel: Dimensionless exclusion volume v ey[ /v e s versus dimension d for 
spherical and cubical hyperplates. Right panel: Lower bounds on the percolation threshold rj c 
versus dimension d for spherical and cubical hyperplates. 

r that is a characteristic length scale of the hyperplate. For spherical hyperplates, r is 
chosen to be radius of the hyperplate. For cubical hyperplates, r is the radius of a (d — 1)- 
sphere whose (d — 1) -dimensional volume (or d- dimensional area) is the same as that of 
the cubical hyperplate. Clearly, v ex /v c g for spherical hyperplates is smaller than that for 
cubical hyperplates in all dimensions d > 3. Lower bounds on r\ c for spherical and cubical 
hyperplates, as obtained from inequality f lT3|) . are shown in the right panel of Fig. HI As 
expected, the lower bound for spherical hyperplates is always greater than that for cubical 
hyperplates in all dimensions d > 3. 



VIII. CONCLUSIONS AND DISCUSSION 



We have exploited the principle that low- dimensional continuum percolation behavior en- 
codes high- dimensional information and recent analytical results on the percolation thresh- 
old 1] C of continuum percolation models [?| to formulate a scaling relation for rj c that is 
applicable to a wide variety of overlapping nonspherical hyperparticles with random orien- 
tations for any dimension d. This scaling relation, which is also an upper bound on rj c , 
depends on the <i-dimensional exclusion volume v ex associated with a hyperparticles and 
was determined analytically for, among other shapes, various polygons for d = 2, Platonic 
solids, sphero cylinders, parallepipeds and zero- volume plates for d = 3 and their appropriate 
generalizations for d > 4. The scaling relation already provides a good estimate of r\ c for 
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FIG. 5: (Color online) Exclusion region associated with equilateral triangles with volume v with 
specific relative particle orientations in M 2 . Left panel: Triangles are perfectly aligned. The 
associated exclusion region is a regular hexagon (yellow region) whose edge length is the same as 
that of the triangle, resulting in v ex /v = 6. Right panel: Triangles are oriented in such a way 
that one can be mapped to another by a center inversion operation. Here two particles always 
form edge-to-edge contacts. The associated exclusion region is a larger equilateral triangle (yellow 
region) whose edge length is twice that of the original triangle, resulting in v ex /v = 4. 

d = 2, becomes increasingly accurate as d increases, and becomes exact in the high-c? limit. 
We also conjectured that overlapping hyperspheres provide the highest threshold among all 
overlapping systems of convex hyperparticle shapes for any d. 

Although the main focus of this paper concerned the percolation threshold of randomly 
oriented overlapping hyperparticles, the bound (j7j) and the scaling relation ( I54~|) are also valid 
for overlapping particles with specific orientational distributions when the general relation 
(Q for the exclusion volume v ex is used. To illustrate this point, we appeal to two-dimensional 
examples depicted in Fig. which shows the exclusion regions associated with equilateral 
triangles in R 2 with two different particle orientations: fully aligned and edge-to-edge . 
Equilateral triangles are not centrally symmetric and so its exclusion volume v ex is minimized 
for the edge-to-edge configuration (i.e., v ex /v = 4) and maximized when the particles are 
aligned (i.e., v ex /v = 6). The exclusion volume for other relative orientations must be 
between these two extremal values. For example, we see from Table EU v ex /v = 5.30797. . . 
for randomly oriented equilateral triangles. Thus, in order of increasing values of v ex /v for 
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these three examples, we have the lower bounds r\ c > 0.25 (edge-to-edge) , r\ c > 0.188395 . . . 
(random), and r\ c > 0.166666... (aligned), respectively, as predicted by (J7J). The actual 
threshold trends for these three cases should follow those given by these bounds, as does the 
corresponding predictions of the scaling relation: r] c « 1.1281, r] c « 0.8501, and r] c « 0.7520. 

It is noteworthy that the size of the dimensionless exclusion volume v ex /v of a convex 
body with some specified orientation distribution is closely related to how densely such 
nonoverlapping objects can be packed in space 33|. Recall that v ex /v for a centrally sym- 
metric convex hyperparticle is minimized when the bodies are aligned Q. For noncentrally 
symmetric convex shapes, nonaligned orientations result in a smaller value of v ex /v relative 
to that of the aligned case (see Fig. [5]). These exclusion- volume properties of centrally and 
noncentrally symmetric bodies play a major role in recently proposed organizing principles 
concerning optimal (i.e., maximally dense) packings of both hard convex and hard concave 



particles in M. 3 33]. Among other results, it was conjectured that the optimal packing of a 



centrally symmetric particle is achieved by the associated optimal Bravais-lattice packing in 
which all the particles are aligned; while the optimal packing of a particle without central 
symmetry is generally given by a non-Bravais-lattice packing in which the particles have 
different nonaligned orientations. 

In the Introduction, we discussed the duality between the equilibrium hard-hyperparticle 
fluid system and the continuum percolation model of overlapping hyperparticles of the same 



shape. A consequence of this duality re 



ation and the so-called decorrelation principle for 



disordered hard-hyperparticle packings [34|, |35| led to the result that the large-c? percola- 



tion threshold 7] c of overlapping hyperparticles is directly related to the large-ci freezing- 
point density (onset of disorder-order phase transition) of corresponding equilibrium hard- 
hyperparticle models Q. In a future work, we shall explore this duality relationship to 
predict the onset of phase transitions in hard-hyperparticle fluids across dimensions. 
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Appendix A: Explicit Calculation of the Mean Width w of A Three-Dimensional 
Spherical Hyperplate in R 4 



In Sec. II, we presented an explicit method to calculate the mean width w = 2R of a 
convex body in M. d . Recall that the width w(n) of a convex body is the smallest distance 
between two impenetrable parallel (d — 1)— dimensional hyperplanes contacting the bound- 
aries of the body, where n is a unit vector orthogonal to the hyperplanes. The mean width 
w is the average of the width w(n) of the body such that n is uniformly distributed over 
the unit sphere S^ 1 G M . Here, we employ this method to obtain w of a three-dimensional 
spherical hyperplate (i.e., 3-sphere) of radius a in R 4 . 

Due to the rotational symmetry of the 3-sphere, its orientation in M 4 is specified by a 
unit vector u, which is analogous to the case of a circular disk in M 3 ; see Table IHIl The 
width w(n) of the 3-sphere is given by 



w 



n) = 2ay/l- < u, n > 2 = 2a sin 0, (A-l) 



where < x, y > is the inner product of two vectors x and y, and 9 is the angle between u 
and n. Therefore, the mean width is given by 

w = — / 2a sin dS 3 = — / / / 2a sin 2 9 sin dOd^dip, (A-2) 
J " Jo Jo Jo 



where 



dS 3 = / / sin sin dOdcpdip = 2vr 2 , (A-3) 
Jo Jo Jo 

is the 4-dimensional solid angle. Carrying out the integration in Eq. flA-2j) yields the mean 
width w of 3-sphere in R 4 , i.e., 

1 167ra 8a 

* = w x — = T, (A - 4) 



Appendix B: Rescaled-Particle Method for Estimating r] c of Randomly Oriented 
Overlapping Platonic Solids 



In Ref. [23|, we devised a highly efficient rescaled-particle method to estimate rj c for d- 
dimensional hyperspheres and oriented hypercubes, utilizing the tightest lower bound on 
7] c , i.e., the pole of the [2, 1] Pade approximant of the mean cluster size S, obtained in 
Ref. |7|; see Eq. (119) in that paper. The basic idea behind this method is to begin 
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with static configurations of overlapping particles at a reduced density t]q that is taken 
to be the best lower-bound value and then to increase rj by rescaling the particle sizes 
until a system-spanning cluster forms. This algorithm enables one to accurately estimate 
r] c in a computationally efficient manner, even in high dimensions 23]. Here, we adapt 
the method to investigate percolation properties of randomly oriented overlapping Platonic 
solids, including tetrahedra, icosahedra, dodecahedra, octahedra, and cubes in M. 3 . Except 
for overlapping cubes, whose percolation properties have been well-studied, the percolation 
thresholds for the other polyhedra have heretofore not been investigated. 

Initially, a Poisson distribution of a large number of points in the periodic simulation 
domain is generated. Each point is then taken to be the centroid of a polyhedron with 
a random orientation and a characteristic length £ (e.g., the diameter of the associated 
circumsphere) . The initial value of £q is chosen so that the initial reduced density r/o of the 
system equals the lower-bound value given by (J7J). For each polyhedron i, a near- neighbor 
list (NNL) is obtained that contains the centroids of the polyhedra j whose distance Dij 
to particle i is smaller than 7^0 (7 > 1)- The value of 7 generally depends on the specific 
polyhedron of interest. A good choice for 7 is a value that yields a NNL list that only 
contains polyhedra that overlap at the threshold. In our simulations, we used 7 G [1.25, 1.5], 



depending on the asphericity of the shapes [36|, |37| . Then the sizes of the polyhedra are slowly 



and uniformly increased by increasing £ Q , leading to an increase of the reduced density by 
an incrementally small amount 5r). After each system rescaling, the polyhedra in the NNL 
are checked for overlap using the separation axis theorem 37[ and the largest cluster is 
identified. This process is repeated until a system-spanning cluster forms. 

Since static particle configurations and a predetermined NNL are used, the complexity of 
cluster identification is significantly reduced. Because initial configurations with relatively 
high reduced density values (i.e., the lower bound values) are used, the amount of rescaling 
before the system percolates is much smaller than that from an initial low-density configu- 
ration. Furthermore, the reduced density rj can be incremented by a small amount 5r) as r) c 
is approached, leading to a more accurate estimate of the percolation threshold. 

For each of the Platonic solids, we use several distinctly different system sizes, i.e., 
N = 5000, 10000, 50000, 100000, and for each system size, we generate 500 indepen- 
dent realizations of the overlapping polyhedra. The results are extrapolated by spline fitting 



35 



the finite-system-size data in a log-log plot to obtain the infinite-system-size estimate of rj c . 
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